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Abstract— We analyze the three-dimensional kinematics of about 82000 Tycho-2 stars 
belonging to the red giant clump (RGC). First, based on all of the currently available 
data, we have determined new, most probable components of the residual rotation vector 
of the optical realization of the ICRS/HIPPARCOS system relative to an inertial frame 
of reference, (w t , u v , oj z ) = (-0.11,0.24,-0.52) ± (0.14,0.10,0.16) mas yr -1 . The stellar 
proper motions in the form fi a cos5 have then be corrected by applying the correction 
uj z = —0.52 mas yr -1 . We show that, apart from their involvement in the general Galactic 
rotation described by the Oort constants A = 15.82 ± 0.21 km s _1 kpc -1 and B = 
—10.87 ± 0.15 km s -1 kpc -1 , the RGC stars have kinematic peculiarities in the Galactic 
yz plane related to the kinematics of the warped stellar-gaseous Galactic disk. We show 
that the parameters of the linear Ogorodnikov-Milne model that describe the kinematics 
of RGC stars in the zx plane do not differ significantly from zero. The situation in the 
yz plane is different. For example, the component of the solid-body rotation vector of the 
local solar neighborhood around the Galactic x axis is Af- = -2.6 ± 0.2 km s -1 kpc -1 . 
Two parameters of the deformation tensor in this plane, namely M+ = 1.0 ± 0.2 km s _1 
kpc -1 and M 33 — M 22 = —1.3 ± 0.4 km s _1 kpc -1 , also differ significantly from zero. On 
the whole, the kinematics of the warped stellar-gaseous Galactic disk in the local solar 
neighborhood can be described as a rotation around the Galactic x axis (close to the line 
of nodes of this structure) with an angular velocity (—3.1 ±0.5) < < (—4.4 ±0.5) km 
s _1 kpc -1 . 



INTRODUCTION 

Analysis of the large-scale structure of neutral hydrogen revealed a warp of the gaseous 
disk in the Galaxy (Westerhout 1957; Burton 1988). The results of studying this structure 
based on currently available data on the HI and HII distributions are presented in Kalberla 
and Dedes (2008) and Cersosimo et al. (2009), respectively. This structure is revealed by 
the spatial distribution of stars and dust (Drimmel and Spergel 2001), by the distribution 
of pulsars in the Galaxy (Yusifov 2004), by HIPPARCOS OB stars (Miyamoto and Zhu 
1998), and by the 2 MASS red giant clump (Momany et al. 2006). 

Several models were suggested to explain the nature of the Galactic disk warp: (1) 
the interaction between the disk and a nonspherical dark matter halo (Sparke and Caser- 
tano 1988); (2) the gravitational influence of the Galaxy's nearest satellites (Bailin 2003); 
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(3) the interaction of the disk with a near-Galaxy flow formed by high-velocity hydro- 
gen clouds that resulted from mass exchange between the Galaxy and the Magellanic 
Clouds (Olano 2004); (4) an intergalactic flow (Lopez-Corredoira et al. 2002); and (5) 
the interaction with the intergalactic magnetic field (Battaner et al. 1990). 

By analyzing nearby stars from HIPPARCOS (1997), Dehnen (1998) showed that the 
distribution of their residual velocities in the V y — V z plane agreed satisfactorily with 
various rotation models of the warped disk. Miyamoto et al. (1993) and Miyamoto and 
Zhu (1998) determined the rotation parameters of the warped stellar-gaseous disk by 
analyzing giant stars of various spectral types and a sample of HIPPARCOS 0-B5 stars. 
Thus, there is positive experience in solving this problem using data on stars relatively 
close to the Sun. 

Studying the three-dimensional kinematics of stars requires that the observational 
data be free from the systematic errors related to the referencing of the optical real- 
ization of the ICRS/HIPPARCOS system to the inertial frame of reference specified by 
extragalactic sources. The modern standard system of astronomical coordinates, ICRS 
(International Celestial Reference System), is realized by the catalog of positions for 212 
compact extragalactic radio sources uniformly distributed over the entire sky observed by 
the radio interferometry technique (Ma et al. 1998). In the optical range, the first real- 
ization of the ICRS was the HIPPARCOS catalog. The application of various methods of 
analysis shows that there is a small residual rotation of the ICRS/HIPPARCOS system 
relative to the inertial frame of reference with tu z —0.4 ± 0.1 mas yr _1 (Bobylev 2004a, 
2004b). At present, there are several new results of comparing individual programs with 
the catalogs of the ICRS/HIPPARCOS system. One of our goals is to determine the 
most probable components of the residual rotation vector of the optical realization of the 
ICRS/HIPPARCOS system relative to the inertial frame of reference. 

The main goal of this paper is to study the local kinematics of the warped stel- 
lar-gaseous Galactic disk by analyzing the motions of Tycho-2 stars that belong to the 
red giant clump (RGC). Occupying a compact region on the Hertzsprung-Russell diagram, 
these stars are a kind of "standard candles" . Their estimated photometric distances are 
known with a mean accuracy of at least 20-30%. RGC stars are distributed uniformly in 
the spatial region and over the celestial sphere, which is an important property when the 
three-dimensional spatial motions of stars are analyzed. 

DATA 

The characteristic clump, commonly called the red giant clump, on the 
Hertzsprung-Russell diagram is formed mostly by giants with masses from 3M Q to 9M . 
These stars spend the bulk of their lifetime on the main sequence as B-type stars. At the 
core helium burning stage, they evolve toward the RGC almost without changing their 
luminosity. At the RGC stage, despite the slow change in color, the luminosity of such a 
star remains almost constant. Therefore, RGC stars are convenient as "standard candles" 
for distance determinations. The RGC also incorporates other giants at different evolu- 
tionary stages. In addition, owing to probabilistic selection methods using reduced proper 
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Fig. 1. Color-absolute magnitude diagram for 82324 RGC stars from the range of dis- 
tances 0.3-1 kpc (a) and the sample of nearby RGC stars with reliable (e^/n < 10%) 
HIPPARCOS parallaxes (b). The isochrones (Yi et al. 2003) for ages of 1, 5, and 10 Gyr 
with a nearly solar metallicity, Z = 0.02 (dashed lines), and a higher metallicity, Z = 0.04 
(solid lines), are plotted. 
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motions, the RGC sample can be diluted by various stars from adjacent regions of the 
Hertzsprung-Russell diagram, both supergiants and dwarfs. According to the estimates 
by various authors, the admixture is small, may be 10-15% (Rybka 2006; Gontcharov 
2008). 

We used the list of 97000 RGC stars selected by Gontcharov (2008) by the dereddened 
color index and the reduced proper motion based on Tycho-2 (Hog et al. 2000) and 2MASS 
(Skrutskie et al. 2006) data. These stars occupy the region with 0. m 5 < J — Ks < m .9 
and — l m < Mk s < — 2 m on the corresponding diagram. For their selection using the 
reduced proper motion, calibration based on HIPPARCOS stars with the most reliable 
data was performed. Details of the procedure are described in Gontcharov (2008). An 
estimate of the photometric distance and interstellar extinction is available for each star. 

As was shown by Bobylev et al. (2009), the kinematic parameters of the Ogorod- 
nikov-Milne model for relatively close (r < 0.2 — 0.3 kpc) and distant (r > 1 kpc) stars 
are determined with large errors when using the photometric distances of RGC stars. 
Therefore, here we use 82 324 stars from the range of distances 0.3 kpc < r < 1 kpc, with 
the mean distance being r = 0.57 ±0.17 kpc. To eliminate the halo stars from the sample, 
the following constraint on the proper motions was applied: ^ (fx a cos 5) 2 + (l^s) 2 < 300 
mas yr _1 . 

As was shown by Bobylev et al. (2009), from a statistical analysis of the velocity 
dispersions, we found that about 20% of the RGC stars are very young, while an over- 
whelming majority of the remaining stars are characterized by the kinematics of thin-disk 
stars. We see from Fig. 1, where the color-absolute magnitude diagram is presented with 
a grid of isochrones with various ages and metallicities, that the RGC stars (Fig. la) 
are fairly young. Figure lb presents about 650 stars from Gontcharov's catalog selected 
as RGC ones with parallaxes from the HIPPARCOS catalog (with a parallax error of 
less than 10%). Although the segment of the isochrones near the RGC is very sensitive 
to metallicity, we clearly see from Fig. lb that, first, the RGC stars proper have ages 
younger than 1 Gyr and, second, the admixture of dwarfs is insignificant. 

THE OGORODNIKOV-MILNE MODEL 

We use a rectangular Galactic coordinate system with the axes directed away from the 
observer toward the Galactic center (/=0°, 6=0°, the X axis or axis 1), along the Galactic 
rotation (Z=90°, 6=0°, the Y axis or axis 2), and toward the North Galactic Pole (6=90°, 
the Z axis or axis 3). 

In the linear Ogorodnikov-Milne model (Ogorodnikov, 1965), we use the notation 
introduced by Clube (1972, 1973) and used, for example, by Vityazev and Tsvetkov 
(2009). 

The observed velocity V(r) of a star with a heliocentric radius vector r, is described, 
to the terms of the first order of smallness t/Rq <C 1, by the equation in vector form 

V(r) = V + Mr + V, (1) 
where V Q (X Q ,Y Q , Z Q ) is the peculiar velocity of the Sun relative to the centroid of the 
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stars under consideration; V is the residual velocity of the star (here, the residual stellar 
velocities are assumed to have a random distribution); M is the displacement matrix 
(tensor) whose components are the partial derivatives of the velocity u(ui, u 2 , w 3 ) with 
respect to the distance r(ri,r 2 ,r 3 ), where u = V(i?) — \(Rq), and R and R are the 
Galactocentric distances of the star and the Sun, respectively. Then, 

M„= (j^j , (p, 9 = 1,2,3). (2) 

All nine elements of the matrix M can be determined using the three components of 
the observed velocities-the stellar radial velocities and proper motions. Having only the 
proper motions, we can write the conditional equations 

4.74r/i^cos6 = X sin/ — F cos/ + (3) 
+r [— cos 6 cos I sin IM U — cos 6 sin 2 IM\ 2 — sin b sin /Mi 3 + cos b cos 2 /M 2] + 
+ cos b sin / cos IM 22 + sin b cos /M 23 ] , 
4. 74rfi b = X Q cos / sin b + F sin / sin b — Z & cos 6 + (4) 
+r [— sin 6 cos b cos 2 /M n — sin b cos 6 sin I cos IM 12 — 
— sin 2 b cos /Mi 3 — sin b cos 6 sin / cos IM 21 — sin 6 cos b sin 2 /M 22 — sin 2 b sin /M 23 + 
+ cos 2 b cos /M 31 + cos 2 b sin /M 32 + sin b cos 6M 33 ] , 

from which it follows that the terms should be grouped in several cases. It is useful to 
divide the matrix M into symmetric, M + (local deformation tensor), and antisymmetric, 
(rotation tensor), parts: 

This allows the conditional equations to be written as 

4.74?t^cos& = X Q sin/ — F Q cos/ + (6) 
+r[—M~ 2 cos / sin b — M~ sin I sin b + M~ cos 6+ 
+M+ cos 21 cos 6 - M+ sin I sin 6 + M+ cos Z sin 6 - 0.5(M+ - M+) sin 2/ cos 6], 

4.74rfi b = X Q cos / sin 6 + F sin I sin 6 — Z cos b + (7) 
+r[M 3 - sin I - M~ cos I - 0.5M+ sin 2/ sin 2b + M+ cos / cos 26+ 
+M+ sin I cos 26 - 0.5(M+ - M+) cos 2 Z sin 2b + 0.5(M+ - M+) sin 26]. 

Equation (6) was obtained from Eq. (3) by adding two terms to its right-hand side, 0.5M 31 
and 0.5M 32 , and subtracting them. 

The system of equations (6) and (7) then becomes very convenient for its simultaneous 
solution. Indeed, as can be seen from Eq. (6), the two pairs of unknowns to be determined, 
M~ and M+, along with M~ and M+, have identical coefficients, sin /sin 6 and cos /sin 6, 
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respectively. In this case, the variables cannot be separated and can be found only from 
the simultaneous solution of the system of equations (6) and (7). 

In addition, we see that one of the diagonal terms of the local deformation tensor 
remains uncertain. Therefore, we determine differences of the form (M+ — M+) and 
(M+-M+). 

The quantities M 32 , M 13 , M 12 , are the components of the solid-body rotation vector of 
a small solar neighborhood around the x, y, z axes, respectively. According to our chosen 
rectangular coordinate system, the rotations from axis 1 to 2, from axis 2 to 3, and from 
axis 3 to 1 are positive. M~ is equivalent to the Oort constant B. Each of the quantities 
M+,M+,M+ describes the deformation in the corresponding plane. M+ is equivalent to 
the Oort constant A. 

The diagonal components of the local deformation tensor M+, M+, M+ describe the 
general local contraction or expansion of the entire stellar system (divergence). The 
system of conditional equations (6) and (7) contain eleven unknowns to be determined 
by the least-squares method. In addition to the simultaneous solution of the system of 
equations (6) and (7), here we analyze the results of their separate solutions. In this case, 
Eq. (7) remains unchanged, while Eq. (6) is reduced to the form 

4.74r/ii cos b = X Q sin / - Y Q cos / + (8) 
+r [M 23 cos I sin b — M 13 sin / sin b + M~ cos b + M+ cos 21 cos b— 
-0.5(M+ - M+) sin 21 cos 6], 

where there are only seven independent variables. 

INERTIALITY OF THE OPTICAL ICRS/HIPPARCOS SYS- 
TEM 

Table 1 presents all of the currently known results of comparing individual programs with 
the catalogs of the ICRS/HIPPARCOS system. 

First, we will briefly describe the results, most of which were used by Kovalevsky et al. 
(1997) to calibrate the HIPPARCOS catalog and to estimate its residual rotation relative 
to the system of extragalactic sources and by Bobylev (2004b) to solve this problem. 

(1) The NPM1 solution. We used the result of comparing the stellar proper motions 
from the NPM1 (Klemola et al. 1994) and HIPPARCOS catalogs performed by the 
Heidelberg group (Kovalevsky et al. 1997). This solution was obtained in the range of 
magnitudes 10 m 5-ll m 5, where (Fig. 1 in Platais et al. 1998a) the HIPPARCOS-NPM1 
stellar proper motion differences have a "horizontal" pattern near zero. In our opinion, 
the NPM1 proper motions in this magnitude range are free from the influence of the 
magnitude equation, which is significant in this catalog, to the greatest extent. 

(2) The NPM2 solution. The stellar proper motions from the NPM2 and HIPPARCOS 
catalogs were compared by Zhu (2003). However, there are no images of galaxies on NPM2 
photographic plates and, hence, these proper motions are relative. Therefore, we do not 
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Table 1: Components of the residual rotation vector of the optical realization of the 
ICRS/HIPPARCOS system relative to the inertial frame of reference 



Method 


Px / Py / Pz 


N* 


N 

area 


u x , mas yr 1 


Uy, mas yr 1 


u z , mas yr 1 


NPM1 


10/5/16 


2616 


899 


-0.76 ±0.25 


±0.17 ±0.20 


-0.85 ±0.20 


NPM2 


(*) 


3519 


347 


-0.11 ±0.20 


-0.19 ±0.20 


-0.75 ±0.28 


SPM2 


22/9/28 


9356 


156 


±0.10 ±0.17 


±0.48 ±0.14 


-0.17 ±0.15 


Kiev 


1/1/1 


415 


154 


-0.27 ±0.80 


±0.15 ±0.60 


-1.07 ±0.80 


Potsdam 


2/1/3 


256 


24 


±0.22 ±0.52 


±0.43 ±0.50 


±0.13 ±0.48 


Bonn 


5/3/6 


88 


13 


±0.16 ±0.34 


-0.32 ±0.25 


±0.17 ±0.33 


HST 


0.1/0.1/0.05 


78 




-1.60 ±2.87 


-1.92 ±1.54 


±2.26 ±3.42 


EOP 


8/2/- 






-0.93 ±0.28 


-0.32 ±0.28 




PUL2 


3/1/4 


1004 


147 


-0.98 ±0.47 


-0.03 ±0.38 


-1.66 ±0.42 


XPM 


28/9/33 


1 x 10 6 


1431 


-0.06 ±0.15 


+0.17 ±0.14 


-0.84 ±0.14 


VLBI-07 


6/1/5 


46 




-0.55 ±0.34 


-0.02 ±0.36 


+0.41 ±0.37 


Minor PI 


25/5/6 


116 




±0.12 ±0.08 


±0.66 ±0.09 


-0.56 ±0.16 


Mean 1 
Mean 2 








-0.22 ±0.19 
-0.11 ±0.14 


±0.14 ±0.10 
±0.24 ±0.10 


-0.49 ± 0.23 
-0.52 ±0.16 



Note. (*) — the NPM2 solution is not used, iV* is the number of stars/asteroids, N area is the 
number of areas on the celestial sphere, mean 1 is a simple mean (without HST), mean 2 is a 
weighted mean. 



use this solution to derive the mean values of u x , oo y , oo z . It is given in Table 1 to emphasize 
its similarity to the NPM1 solution. 

(3) The SPM2 solution. The stellar proper motions from the SPM2 (Platais et al. 
1998b) and HIPPARCOS catalogs were compared by Zhu (2001). 

(4) The PUL2 solution. The parameters u x , u y , u z were found by comparing the 
Pulkovo PUL2 photographic catalog (Bobylev et al. 2004) and HIPPARCOS. 

(5) The KIEV solution. The stellar proper motions from the GPM1 (Rybka and 
Yatsenko 1997) and HIPPARCOS catalogs were compared by Kislyuk et al. (1997). 

(6) The POTSDAM solution. The parameters u x , u> y , u z of the Potsdam program were 
taken from Hirte et al. (1996). 

(7) The BONN solution. The results of the Bonn program are presented in Geffert et 
al. (1997) and Tucholke et al. (1997). 

(8) The EOP solution. The results of the analysis of Earth orientation parameters 
(EOP) were taken from Vondrak et al. (1997). Only two rotation parameters, u x , and uj y 
are determined in this method. 

(9) The HST solution. The results of stellar observations with the Hubble Space 
Telescope (HST) were taken from Hemenway et al. (1997). 

Now, we will point out several new results that have not been used previously to solve 
this problem. 
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(10) The XPM solution. The XPM catalog (Fedorov et al, 2009) contains absolute 
proper motions for about 275 million stars fainter than 12m derived by comparing their 
positions in the 2MASS and USNO-A2.0 (Monet 1998) catalogs. The absolutization 
was made using about 1.5 million galaxies from the 2MASS catalog of extended sources. 
Thus, the XPM catalog is an independent realization of the inertial frame of reference. 
The stellar proper motions from the XPM and UCAC2 (Zacharias et al. 2004) catalogs 
were compared by Bobylev et al. (2010). The parameters u x ,u y ,u z were calculated using 
about 1 million stars. Among all of the programs listed in Table 1, the XPM solution is 
unique in that it was obtained from the differences of stars covering the entire sky almost 
completely, except for the zone 5 > 54°, where there are no UCAC2 stars. 

(11) The VLBI-07 solution. Boboltz et al. (2007) analyzed the positions and proper 
motions of 46 radio stars and obtained new parameters of the mutual orientation of the 
optical realization (HIPPARCOS) and the radio system. 

(12) The "MINOR PLANETS" solution. Chernetenko (2008) estimated the rotation 
parameters of the HIPPARCOS system relative to the DE403 and DE405 coordinate 
systems of ephemerides by analyzing a long-term series of asteroid observations. This 
result suggests that either the dynamical DE403 and DE405 theories need to be improved 
or the HIPPARCOS system needs to be corrected. We reduced the weight of this solution 
by half because of the possible contribution from the inaccuracy of the dynamical DE403 
and DE405 theories. 

The weight of each of the comparison catalogs was taken to be inversely proportional 
to the square of the random error e w in the corresponding quantities u x ,u y ,u z and was 
calculated from the formula 



Not all of the authors use the equations to determine u! x ,ui y , oo z in the form in which they 
were suggested by Lindegren and Kovalevsky (1995): 



where the catalog-HIPPARCOS differences are on the left-hand sides of the equations. 
Therefore, in several cases, the signs of the quoted quantities were reduced to the necessary 
uniform form (Zhu 2001, 2003; Boboltz et al. 2007). 

The last rows of Table 1 give Mean 1 calculated as a simple mean and Mean 2 that 
was calculated as a weighted mean and is the main result of our analysis. 

Denote the components of the rotation vector around the rectangular equatorial axes 
by u x ,u y ,u z ; then, 




(9) 



Afi a cos 5 = u x cos a sin 5 + u y sin a sin 5 — u z cos S, 
Afis = —u x sin a + Uy cos a, 



(10) 
(11) 



( 




a, 



) 



G 



( 



M, 
M] 



21 



32 



13 



) 




(12) 
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where 



G = 



/ -0.0548 
-0.8734 
V -0.4838 



+0.4941 
-0.4448 
+0.7470 



-0.8677* 

-0.1981 

+0.4560 



(13) 



is the well-known transformation matrix between the unit vectors of the Galactic and 
equatorial coordinate systems. From Eqs. (9) and (10), it is easy to see the relationship 

Assuming the components in the mean-2 solution to be zero, 



between M 13 and tu z 




4.74 



I 


\0.52 



(14) 



for the case where the left- and right-hand sides are expressed in km s kpc 



KINEMATICS OF THE DISK WARP 

The results of determining the kinematic parameters of the Ogorodnikov-Milne model 
using a sample of RGC stars derived by simultaneously solving the system of equations 
(6) and (7) are presented in Table 2. The second column gives the solution obtained 
without applying any corrections to the input data; the third column gives the solution 
for the case where the stellar proper motions in the form fx a cos 5 were corrected by 
applying the correction oj z = —0.52 mas yr™ 1 (14) using Eq. (10). 

The parameters derived by separately solving Eqs. (8) and (7) are presented in Table 

3. 

As can be seen from Table 2, applying the correction affected only three components 
of the rotation tensor: insignificantly M~, noticeably M 3 ~, and most strongly M~ 3 , which 
is explained by the structure of the matrix G (13). Since the components M~ and M+ 
found (the third column of Table 2) do not differ significantly from zero, they may be set 
equal to zero. 

Suppose that the values of M 2 ~ and M+ (the Oort constants) found describe only the 
rotation around the Galactic z axis, while the motion in the yz plane is independent. 
Let us now consider the displacement tensor M w that describes the kinematics in the yz 
plane: 



w 



M 22 
M 32 



M 23 
M 33 



/ du 2 du 2 \ 
~dr2~ ~dr^ 
dus <9«3 



(15) 



As has already been noted, when using only the stellar proper motions, we can determine 
only the difference (M+ - M+). The identity (M+ - M+) = (M 33 - M 22 ) is valid for 
the diagonal elements. As can be seen from Table 2, the value of this quantity differs 
significantly from zero. 

Three cases are possible when analyzing tensor (15): 
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1) M 22 ^ 0, M 33 = 0; 

2) M 22 = 0, M 33 ^ 0; 

3) M 22 = 0, M 33 = 0. 

For the completeness of the picture, note that case 4 is also possible: M 22 ^ 0, M 33 7^ 0; 
since this requires data on the stellar radial velocities, this case is not considered here. 

Consider case 1, M 33 = 0,, using the data from the third column of Table 2, M 22 = 
1.3 ± 0.4 km s _1 kpc -1 . The components of the displacement tensor Mw, the symmetric 
deformation tensor Myy, and the antisymmetric rotation tensor are (km s _1 kpc -1 ) 

The deformation tensor in the principal axes is (km s" 1 kpc -1 ): 



Ai 
A 2 
1.8 
-0.5 

and the angle between the positive direction of the Oy axis and the first principal axis of 
this ellipse is 29 ± 6°. 

Equation (1) can be written as 

V = V + grad F + (wxr). 

Analysis of grad F shows (Sedov 1970) that an infinitesimal sphere of radius r composed 
of points of the medium at time t, 

2,2,2 2 
x + y + z = r 

transforms into a deformation ellipsoid after time At, 

x* 2 y* 2 z* 2 9 



(1 + AxAt) 2 (l + A 2 At) 2 (1 + A 3 A*) : 

which is shown in Fig. 2b for our two-dimensional case. 

For case 2, we have M 33 = —1.3 ± 0.4 km s -1 kpc" 1 and M 22 = 0. The deformation 
tensor in the principal axes is (km s -1 kpc -1 ) 



M+ = 



0.5 
-U 
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Table 2: Kinematic parameters of the Ogorodnikov-Milne model found by simultaneously solv- 
ing the system of equations (6) and (7) 



Parameter 


Without correction 


With correction 


x & 


7.99 ±0.10 


7.99 ±0.10 




16.40 ±0.10 


16.40 ±0.10 


z Q 


6.72 ±0.09 


6.72 ±0.09 


A = M+ 


15.82 ±0.21 


15.82 ±0.21 




-1.40 ±0.18 


-2.60 ±0.18 


Mf 3 


-1.99 ±0.18 


-0.15 ±0.18 


B = M 2 1 


-11.99 ±0.15 


-10.87 ±0.15 


Mtl-22 


-7.75 ±0.39 


-7.75 ±0.39 


M+ 


-0.47 ±0.23 


-0.47 ±0.23 




1.00 ±0.22 


1.00 ±0.22 


-^33-22 


-1.26 ±0.44 


-1.26 ±0.44 



Note: Xq,Yq, Z@ in km s , other parameters in km s kpc 



Table 3: Kinematic parameters of the Ogorodnikov-Milne model found by separately solving 
Eqs. (8) and (7) 



Parameter 


Without correction" 


With correction 61 


Without correction b 


With correction 6 


X & 


7.89 ±0.12 


7.89 ±0.12 


8.51 ±0.21 


8.51 ±0.21 


Y Q 


16.13 ±0.13 


16.13 ±0.13 


17.40 ± 0.21 


17.40 ± 0.21 


Z & 






6.73 ± 0.08 


6.73 ± 0.08 


A = M+ 


15.76 ±0.25 


15.76 ±0.25 


15.99 ± 0.48 


15.99 ±0.48 


^32 






-1.47 ±0.25 


-2.66 ±0.25 


M 13 






-1.46 ±0.25 


0.38 ±0.25 


B = 


-12.02 ±0.17 


-10.91 ±0.17 






Kl-22 


-8.21 ±0.46 


-8.21 ±0.46 


-4.31 ± 1.03 


-4.31 ± 1.03 


M+ 






0.17 ±0.34 


0.17 ±0.34 


+ ^2 + 3 






1.11 ±0.32 


1.11 ±0.32 


-^33-22 






0.29 ± 0.59 


0.29 ±0.59 


M 23 


2.24 ±0.46 


3.43 ± 0.46 






M13 


-2.93 ± 0.43 


-1.09 ±0.43 







a The parameters were found by solving Eq. (8). 
fe The parameters were found by solving Eq. (7). 
Note: Xq, Yq, Zq in km s _1 , other parameters in km s _1 kpc 
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the angle between the positive direction of the Oy axis and the first principal axis of this 
ellipse is 29 ± 6°. 

For case 3, where M 33 = and M 22 = 0, the deformation tensor has two roots: 
Ai = 1.0 km s -1 kpc -1 and A2 = —1.0 km s _1 kpc -1 ; therefore, the first principal axis of 
this ellipse is oriented at an angle of 45° to the Oy axis. 

In all three cases, the divergence 0.5(M 22 + M 33 ) is insignificant. It is +0.6 ± 0.4, and 
km s" 1 kpc -1 for the first, second, and third cases, respectively. 

As a result, we can conclude that the kinematics in the yz plane can be described as 
a rotation around the Galactic x axis with an angular velocity flw = M~ 2 — X±. This 
angular velocity is —4.4 ± 0.5, —3.1 ± 0.5, and —3.6 ± 0.3 km s -1 kpc -1 for the first, 
second, and third cases, respectively. 

Let us estimate the linear velocities in the yz plane. The mean photometric distance 
for our sample of RGC stars is 7 = 0.57 ± 0.17 kpc. The linear solidbody rotation 
velocity is then M~ 2 ■ r = —1.5 ± 0.1 km s -1 . The maximum deformation velocity is 
Ai • r — 1.0 ± 0.3 km s" 1 and the highest (among the cases considered) linear velocity is 
VL w -r = -2.5 ±0.3 km s _1 . 

Figure 2 gives a qualitative picture that reflects the pattern of our results. The dis- 
tribution of rotation velocity vectors in the yz plane for the case of rotation we found 
is shown in Fig. 2a. Figure 2b shows how a circumference of unit radius turns into a 
deformation ellipse after a unit time interval. 

According to the results of our separate solution of Eqs. (8) and (7) (the third and 
fifth columns of Table 3), for example, for case 1, we have 



We can see from our comparison of matrixes (19) and (16), (20) and (17), (21) and (18) 
that there are differences only in the diagonal element. This leads us to conclude that both 
approaches used (the simultaneous and separate solution of the equations) are equivalent, 
because they yield coincident results. 

When writing Eq. (14), we set the components equal to zero, u x = and oo y = 0. As 
we see from the Mean 2 solution in Table 1, the value of the component u y was found 
to be nonzero outside the 2a range. However, the kinematic parameters calculated by 
taking this into account, i.e., u y 7^ and u z 7^ 0, have no significant differences between 
the solutions reflected in Tables 2 and 3 and, hence, we do not give them. 



Drimmel et al. (2000) provide arguments for the model of precession of the warped disk in 
the zx plane (rotation around the y axis) with an angular velocity of —25 km s -1 kpc -1 . 




(21) 



(19) 



(20) 



DISCUSSION 
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However, as can be seen from Table 2, as a result of applying the correction oo z = —0.52 
mas yr _1 , the component of the rotation tensor around the y axis (M~) has an almost 
zero value and the component of the deformation tensor in this plane (M+) does not differ 
significantly from zero. 

The angular velocity Qw estimated here is inconsistent with the analysis of the motion 
of giants of various spectral types (O-M) performed by Miyamoto et al. (1993), who found 
a positive direction of rotation around the x axis (directed toward the Galactic center). 
Note that the analyzed stellar proper motions were determined in the FK5 (Fundamental 
Catalog 5) system, which is noticeably distorted by the uncertainty in the precession 
constant. On the other hand, having analyzed the proper motions of 0-B5 stars in 
the HIPPARCOS system, Miyamoto and Zhu (1998) also reached the conclusion about a 
positive rotation around the 

At present, it is impossible to choose between the three forms of the deformation tensor 
considered in the previous section due to the absence of data on stellar radial velocities 
and distances. Note that based on a sample of 3632 relatively close RGC stars with known 
HIPPARCOS trigonometric parallaxes and radial velocities, we estimated the parameters 
M 22 = -1.3±3.2 km s" 1 kpc" 1 and M 33 = -1.1 ±2.6 km s" 1 kpc" 1 (Bobylev et al. 2009). 
Consequently, the difference (M+ — M+) is very close to zero and, hence, the orientation 
of the deformation ellipse is close to 45°, as in case 3 considered. 

In the Introduction, we listed various hypotheses about the nature of the Galactic 
disk warp. In our opinion, the hypothesis about an outflow of gas from the Magellanic 
Clouds (Olano 2004) has been worked out in greatest detail. Assuming the distance to 
the Large Magellanic Cloud to be t^mc — 50 kpc, we find that the linear velocity of the 
rotation we revealed at this distance, (—3.1 ± 0.5) < Qw < (—4.4 ± 0.5) km s _1 kpc -1 , 
is \Qw\ ■ f'LMC = (155 4- 220) ± 25 km s _1 . Such velocities are in agreement with the 
estimate of the motion of high- velocity hydrogen clouds relative to the local standard of 
rest, ~ 200 km s" 1 (Olano 2004), and the rotation direction we found is in good agreement 
with the direction of motion of the Magellanic Clouds around the Galaxy during the past 
570 Myr (Fig. 6 from Olano 2004). 

CONCLUSIONS 

We analyzed the kinematics of about 82000 RGC stars from the Tycho-2 catalog. These 
stars are a variety of "standard candles", because they occupy a compact region on the 
Hertzsprung-Russell diagram. Therefore, photometric distance estimates with an accu- 
racy of at least 30% are available for them. The RGC stars considered lie in the range of 
heliocentric distances 0.3 kpc < d < 1 kpc. 

Since these stars are distributed fairly uniformly in space, they are of great interest 
not only in studying the Galaxy's general rotation described by the rotation parameters 
in the Galactic xy plane but also in the other two planes, namely yz and zx. 

For a reliable description of their kinematic peculiarities, we should be confident that 
the observational data (Tycho-2 proper motions) are free from the systematic errors re- 
lated to the referencing of the optical realization of the ICRS/HIPPARCOS system to the 
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inertial frame of reference specified by extragalactic sources. 

Here, we gave much attention to this problem. Based on all of the currently available 
data, we determined new, most probable components of the residual rotation vector of 
the optical realization of the ICRS/HIPPARCOS system relative to the inertial frame of 
reference, (u x ,u y ,u z ) = (-0.11,0.24,-0.52) ± (0.14,0.10,0.16) mas yr -1 . This led us to 
conclude that a small correction, u z = —0.52 mas yr -1 , should be applied to the Tycho-2 
stellar proper motions of the form /j, a cos 5. 

By applying this correction, we showed that, apart from their involvement in the 
general Galactic rotation described by the Oort constants A = 15.8 ± 0.2 km s -1 kpc -1 
and B = —10.9 ± 0.2 km s -1 kpc -1 , the RGC stars considered have peculiarities in the 
yz plane related, in our opinion, to the kinematics of the warped Galactic stellar-gaseous 
disk. 

The component of the solid-body rotation of the local solar neighborhood around the 
Galactic x axis is M~ 2 = —2.6 ± 0.2 km s _1 kpc -1 . The parameter of the deformation 
tensor in this plane M+ = 1.0±0.2 km s -1 kpc -1 and the difference M 33 -M 22 = -1.3±0.4 
km s -1 kpc -1 differ significantly from zero. 

On the whole, the kinematics of the warped Galactic stellar-gaseous disk in the local 
solar neighborhood can be described as a rotation around the Galactic x axis directed 
from the Sun toward the Galactic center with an angular velocity (—3.1 ± 0.5) < Qw < 
(—4.4 ±0.5) km s -1 kpc -1 . Thus, the rotation is around an axis close to the line of nodes 
of this structure. In the case where the difference M 33 — M 22 may be set equal to zero, 
the angular velocity Vt w is —3.6 ± 0.3 km s -1 kpc -1 . If M 22 = and M 33 ^ 0, then 
VL W = -3.1 ± 0.5 km s -1 kpc -1 . If M 22 ^ and M 33 = 0, then Q w reaches -4.4 ± 0.5 km 
s -1 kpc -1 . 
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Parameters of the Local Warp of the Stellar- Gaseous 
Galactic Disk from the Kinematics of Tycho-2 
Nearby Red Giant Clump Stars 

V.V. Bobylev 

Pulkovo Astronomical Observatory, Russian Academy of Sciences, St-Petersburg 

Abstract— We analyze the three-dimensional kinematics of about 82000 Tycho-2 stars 
belonging to the red giant clump (RGC). First, based on all of the currently available 
data, we have determined new, most probable components of the residual rotation vector 
of the optical realization of the ICRS/HIPPARCOS system relative to an inertial frame 
of reference, (w t , u v , oj z ) = (-0.11,0.24,-0.52) ± (0.14,0.10,0.16) mas yr -1 . The stellar 
proper motions in the form fi a cos5 have then be corrected by applying the correction 
uj z = —0.52 mas yr -1 . We show that, apart from their involvement in the general Galactic 
rotation described by the Oort constants A = 15.82 ± 0.21 km s _1 kpc -1 and B = 
—10.87 ± 0.15 km s -1 kpc -1 , the RGC stars have kinematic peculiarities in the Galactic 
yz plane related to the kinematics of the warped stellar-gaseous Galactic disk. We show 
that the parameters of the linear Ogorodnikov-Milne model that describe the kinematics 
of RGC stars in the zx plane do not differ significantly from zero. The situation in the 
yz plane is different. For example, the component of the solid-body rotation vector of the 
local solar neighborhood around the Galactic x axis is Af- = -2.6 ± 0.2 km s -1 kpc -1 . 
Two parameters of the deformation tensor in this plane, namely M+ = 1.0 ± 0.2 km s _1 
kpc -1 and M 33 — M 22 = —1.3 ± 0.4 km s _1 kpc -1 , also differ significantly from zero. On 
the whole, the kinematics of the warped stellar-gaseous Galactic disk in the local solar 
neighborhood can be described as a rotation around the Galactic x axis (close to the line 
of nodes of this structure) with an angular velocity (—3.1 ±0.5) < < (—4.4 ±0.5) km 
s _1 kpc -1 . 



INTRODUCTION 

Analysis of the large-scale structure of neutral hydrogen revealed a warp of the gaseous 
disk in the Galaxy (Westerhout 1957; Burton 1988). The results of studying this structure 
based on currently available data on the HI and HII distributions are presented in Kalberla 
and Dedes (2008) and Cersosimo et al. (2009), respectively. This structure is revealed by 
the spatial distribution of stars and dust (Drimmel and Spergel 2001), by the distribution 
of pulsars in the Galaxy (Yusifov 2004), by HIPPARCOS OB stars (Miyamoto and Zhu 
1998), and by the 2 MASS red giant clump (Momany et al. 2006). 

Several models were suggested to explain the nature of the Galactic disk warp: (1) 
the interaction between the disk and a nonspherical dark matter halo (Sparke and Caser- 
tano 1988); (2) the gravitational influence of the Galaxy's nearest satellites (Bailin 2003); 
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(3) the interaction of the disk with a near-Galaxy flow formed by high-velocity hydro- 
gen clouds that resulted from mass exchange between the Galaxy and the Magellanic 
Clouds (Olano 2004); (4) an intergalactic flow (Lopez-Corredoira et al. 2002); and (5) 
the interaction with the intergalactic magnetic field (Battaner et al. 1990). 

By analyzing nearby stars from HIPPARCOS (1997), Dehnen (1998) showed that the 
distribution of their residual velocities in the V y — V z plane agreed satisfactorily with 
various rotation models of the warped disk. Miyamoto et al. (1993) and Miyamoto and 
Zhu (1998) determined the rotation parameters of the warped stellar-gaseous disk by 
analyzing giant stars of various spectral types and a sample of HIPPARCOS 0-B5 stars. 
Thus, there is positive experience in solving this problem using data on stars relatively 
close to the Sun. 

Studying the three-dimensional kinematics of stars requires that the observational 
data be free from the systematic errors related to the referencing of the optical real- 
ization of the ICRS/HIPPARCOS system to the inertial frame of reference specified by 
extragalactic sources. The modern standard system of astronomical coordinates, ICRS 
(International Celestial Reference System), is realized by the catalog of positions for 212 
compact extragalactic radio sources uniformly distributed over the entire sky observed by 
the radio interferometry technique (Ma et al. 1998). In the optical range, the first real- 
ization of the ICRS was the HIPPARCOS catalog. The application of various methods of 
analysis shows that there is a small residual rotation of the ICRS/HIPPARCOS system 
relative to the inertial frame of reference with tu z —0.4 ± 0.1 mas yr _1 (Bobylev 2004a, 
2004b). At present, there are several new results of comparing individual programs with 
the catalogs of the ICRS/HIPPARCOS system. One of our goals is to determine the 
most probable components of the residual rotation vector of the optical realization of the 
ICRS/HIPPARCOS system relative to the inertial frame of reference. 

The main goal of this paper is to study the local kinematics of the warped stel- 
lar-gaseous Galactic disk by analyzing the motions of Tycho-2 stars that belong to the 
red giant clump (RGC). Occupying a compact region on the Hertzsprung-Russell diagram, 
these stars are a kind of "standard candles" . Their estimated photometric distances are 
known with a mean accuracy of at least 20-30%. RGC stars are distributed uniformly in 
the spatial region and over the celestial sphere, which is an important property when the 
three-dimensional spatial motions of stars are analyzed. 

DATA 

The characteristic clump, commonly called the red giant clump, on the 
Hertzsprung-Russell diagram is formed mostly by giants with masses from 3M Q to 9M . 
These stars spend the bulk of their lifetime on the main sequence as B-type stars. At the 
core helium burning stage, they evolve toward the RGC almost without changing their 
luminosity. At the RGC stage, despite the slow change in color, the luminosity of such a 
star remains almost constant. Therefore, RGC stars are convenient as "standard candles" 
for distance determinations. The RGC also incorporates other giants at different evolu- 
tionary stages. In addition, owing to probabilistic selection methods using reduced proper 
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Fig. 1. Color-absolute magnitude diagram for 82324 RGC stars from the range of dis- 
tances 0.3-1 kpc (a) and the sample of nearby RGC stars with reliable (e^/n < 10%) 
HIPPARCOS parallaxes (b). The isochrones (Yi et al. 2003) for ages of 1, 5, and 10 Gyr 
with a nearly solar metallicity, Z = 0.02 (dashed lines), and a higher metallicity, Z = 0.04 
(solid lines), are plotted. 
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motions, the RGC sample can be diluted by various stars from adjacent regions of the 
Hertzsprung-Russell diagram, both supergiants and dwarfs. According to the estimates 
by various authors, the admixture is small, may be 10-15% (Rybka 2006; Gontcharov 
2008). 

We used the list of 97000 RGC stars selected by Gontcharov (2008) by the dereddened 
color index and the reduced proper motion based on Tycho-2 (Hog et al. 2000) and 2MASS 
(Skrutskie et al. 2006) data. These stars occupy the region with 0. m 5 < J — Ks < m .9 
and — l m < Mk s < — 2 m on the corresponding diagram. For their selection using the 
reduced proper motion, calibration based on HIPPARCOS stars with the most reliable 
data was performed. Details of the procedure are described in Gontcharov (2008). An 
estimate of the photometric distance and interstellar extinction is available for each star. 

As was shown by Bobylev et al. (2009), the kinematic parameters of the Ogorod- 
nikov-Milne model for relatively close (r < 0.2 — 0.3 kpc) and distant (r > 1 kpc) stars 
are determined with large errors when using the photometric distances of RGC stars. 
Therefore, here we use 82 324 stars from the range of distances 0.3 kpc < r < 1 kpc, with 
the mean distance being r = 0.57 ±0.17 kpc. To eliminate the halo stars from the sample, 
the following constraint on the proper motions was applied: ^ (fx a cos 5) 2 + (l^s) 2 < 300 
mas yr _1 . 

As was shown by Bobylev et al. (2009), from a statistical analysis of the velocity 
dispersions, we found that about 20% of the RGC stars are very young, while an over- 
whelming majority of the remaining stars are characterized by the kinematics of thin-disk 
stars. We see from Fig. 1, where the color-absolute magnitude diagram is presented with 
a grid of isochrones with various ages and metallicities, that the RGC stars (Fig. la) 
are fairly young. Figure lb presents about 650 stars from Gontcharov's catalog selected 
as RGC ones with parallaxes from the HIPPARCOS catalog (with a parallax error of 
less than 10%). Although the segment of the isochrones near the RGC is very sensitive 
to metallicity, we clearly see from Fig. lb that, first, the RGC stars proper have ages 
younger than 1 Gyr and, second, the admixture of dwarfs is insignificant. 

THE OGORODNIKOV-MILNE MODEL 

We use a rectangular Galactic coordinate system with the axes directed away from the 
observer toward the Galactic center (/=0°, 6=0°, the X axis or axis 1), along the Galactic 
rotation (Z=90°, 6=0°, the Y axis or axis 2), and toward the North Galactic Pole (6=90°, 
the Z axis or axis 3). 

In the linear Ogorodnikov-Milne model (Ogorodnikov, 1965), we use the notation 
introduced by Clube (1972, 1973) and used, for example, by Vityazev and Tsvetkov 
(2009). 

The observed velocity V(r) of a star with a heliocentric radius vector r, is described, 
to the terms of the first order of smallness t/Rq <C 1, by the equation in vector form 

V(r) = V + Mr + V, (1) 
where V Q (X Q ,Y Q , Z Q ) is the peculiar velocity of the Sun relative to the centroid of the 
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stars under consideration; V is the residual velocity of the star (here, the residual stellar 
velocities are assumed to have a random distribution); M is the displacement matrix 
(tensor) whose components are the partial derivatives of the velocity u(ui, u 2 , w 3 ) with 
respect to the distance r(ri,r 2 ,r 3 ), where u = V(i?) — \(Rq), and R and R are the 
Galactocentric distances of the star and the Sun, respectively. Then, 

M„= (j^j , (p, 9 = 1,2,3). (2) 

All nine elements of the matrix M can be determined using the three components of 
the observed velocities-the stellar radial velocities and proper motions. Having only the 
proper motions, we can write the conditional equations 

4.74r/i^cos6 = X sin/ — F cos/ + (3) 
+r [— cos 6 cos I sin IM U — cos 6 sin 2 IM\ 2 — sin b sin /Mi 3 + cos b cos 2 /M 2] + 
+ cos b sin / cos IM 22 + sin b cos /M 23 ] , 
4. 74rfi b = X Q cos / sin b + F sin / sin b — Z & cos 6 + (4) 
+r [— sin 6 cos b cos 2 /M n — sin b cos 6 sin I cos IM 12 — 
— sin 2 b cos /Mi 3 — sin b cos 6 sin / cos IM 21 — sin 6 cos b sin 2 /M 22 — sin 2 b sin /M 23 + 
+ cos 2 b cos /M 31 + cos 2 b sin /M 32 + sin b cos 6M 33 ] , 

from which it follows that the terms should be grouped in several cases. It is useful to 
divide the matrix M into symmetric, M + (local deformation tensor), and antisymmetric, 
(rotation tensor), parts: 

This allows the conditional equations to be written as 

4.74?t^cos& = X Q sin/ — F Q cos/ + (6) 
+r[—M~ 2 cos / sin b — M~ sin I sin b + M~ cos 6+ 
+M+ cos 21 cos 6 - M+ sin I sin 6 + M+ cos Z sin 6 - 0.5(M+ - M+) sin 2/ cos 6], 

4.74rfi b = X Q cos / sin 6 + F sin I sin 6 — Z cos b + (7) 
+r[M 3 - sin I - M~ cos I - 0.5M+ sin 2/ sin 2b + M+ cos / cos 26+ 
+M+ sin I cos 26 - 0.5(M+ - M+) cos 2 Z sin 2b + 0.5(M+ - M+) sin 26]. 

Equation (6) was obtained from Eq. (3) by adding two terms to its right-hand side, 0.5M 31 
and 0.5M 32 , and subtracting them. 

The system of equations (6) and (7) then becomes very convenient for its simultaneous 
solution. Indeed, as can be seen from Eq. (6), the two pairs of unknowns to be determined, 
M~ and M+, along with M~ and M+, have identical coefficients, sin /sin 6 and cos /sin 6, 
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respectively. In this case, the variables cannot be separated and can be found only from 
the simultaneous solution of the system of equations (6) and (7). 

In addition, we see that one of the diagonal terms of the local deformation tensor 
remains uncertain. Therefore, we determine differences of the form (M+ — M+) and 
(M+-M+). 

The quantities M 32 , M 13 , M 12 , are the components of the solid-body rotation vector of 
a small solar neighborhood around the x, y, z axes, respectively. According to our chosen 
rectangular coordinate system, the rotations from axis 1 to 2, from axis 2 to 3, and from 
axis 3 to 1 are positive. M~ is equivalent to the Oort constant B. Each of the quantities 
M+,M+,M+ describes the deformation in the corresponding plane. M+ is equivalent to 
the Oort constant A. 

The diagonal components of the local deformation tensor M+, M+, M+ describe the 
general local contraction or expansion of the entire stellar system (divergence). The 
system of conditional equations (6) and (7) contain eleven unknowns to be determined 
by the least-squares method. In addition to the simultaneous solution of the system of 
equations (6) and (7), here we analyze the results of their separate solutions. In this case, 
Eq. (7) remains unchanged, while Eq. (6) is reduced to the form 

4.74r/ii cos b = X Q sin / - Y Q cos / + (8) 
+r [M 23 cos I sin b — M 13 sin / sin b + M~ cos b + M+ cos 21 cos b— 
-0.5(M+ - M+) sin 21 cos 6], 

where there are only seven independent variables. 

INERTIALITY OF THE OPTICAL ICRS/HIPPARCOS SYS- 
TEM 

Table 1 presents all of the currently known results of comparing individual programs with 
the catalogs of the ICRS/HIPPARCOS system. 

First, we will briefly describe the results, most of which were used by Kovalevsky et al. 
(1997) to calibrate the HIPPARCOS catalog and to estimate its residual rotation relative 
to the system of extragalactic sources and by Bobylev (2004b) to solve this problem. 

(1) The NPM1 solution. We used the result of comparing the stellar proper motions 
from the NPM1 (Klemola et al. 1994) and HIPPARCOS catalogs performed by the 
Heidelberg group (Kovalevsky et al. 1997). This solution was obtained in the range of 
magnitudes 10 m 5-ll m 5, where (Fig. 1 in Platais et al. 1998a) the HIPPARCOS-NPM1 
stellar proper motion differences have a "horizontal" pattern near zero. In our opinion, 
the NPM1 proper motions in this magnitude range are free from the influence of the 
magnitude equation, which is significant in this catalog, to the greatest extent. 

(2) The NPM2 solution. The stellar proper motions from the NPM2 and HIPPARCOS 
catalogs were compared by Zhu (2003). However, there are no images of galaxies on NPM2 
photographic plates and, hence, these proper motions are relative. Therefore, we do not 
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Table 1: Components of the residual rotation vector of the optical realization of the 
ICRS/HIPPARCOS system relative to the inertial frame of reference 



Method 


Px / Py / Pz 


N* 


N 

area 


u x , mas yr 1 


Uy, mas yr 1 


u z , mas yr 1 


NPM1 


10/5/16 


2616 


899 


-0.76 ±0.25 


±0.17 ±0.20 


-0.85 ±0.20 


NPM2 


(*) 


3519 


347 


-0.11 ±0.20 


-0.19 ±0.20 


-0.75 ±0.28 


SPM2 


22/9/28 


9356 


156 


±0.10 ±0.17 


±0.48 ±0.14 


-0.17 ±0.15 


Kiev 


1/1/1 


415 


154 


-0.27 ±0.80 


±0.15 ±0.60 


-1.07 ±0.80 


Potsdam 


2/1/3 


256 


24 


±0.22 ±0.52 


±0.43 ±0.50 


±0.13 ±0.48 


Bonn 


5/3/6 


88 


13 


±0.16 ±0.34 


-0.32 ±0.25 


±0.17 ±0.33 


HST 


0.1/0.1/0.05 


78 




-1.60 ±2.87 


-1.92 ±1.54 


±2.26 ±3.42 


EOP 


8/2/- 






-0.93 ±0.28 


-0.32 ±0.28 




PUL2 


3/1/4 


1004 


147 


-0.98 ±0.47 


-0.03 ±0.38 


-1.66 ±0.42 


XPM 


28/9/33 


1 x 10 6 


1431 


-0.06 ±0.15 


+0.17 ±0.14 


-0.84 ±0.14 


VLBI-07 


6/1/5 


46 




-0.55 ±0.34 


-0.02 ±0.36 


+0.41 ±0.37 


Minor PI 


25/5/6 


116 




±0.12 ±0.08 


±0.66 ±0.09 


-0.56 ±0.16 


Mean 1 
Mean 2 








-0.22 ±0.19 
-0.11 ±0.14 


±0.14 ±0.10 
±0.24 ±0.10 


-0.49 ± 0.23 
-0.52 ±0.16 



Note. (*) — the NPM2 solution is not used, iV* is the number of stars/asteroids, N area is the 
number of areas on the celestial sphere, mean 1 is a simple mean (without HST), mean 2 is a 
weighted mean. 



use this solution to derive the mean values of u x , oo y , oo z . It is given in Table 1 to emphasize 
its similarity to the NPM1 solution. 

(3) The SPM2 solution. The stellar proper motions from the SPM2 (Platais et al. 
1998b) and HIPPARCOS catalogs were compared by Zhu (2001). 

(4) The PUL2 solution. The parameters u x , u y , u z were found by comparing the 
Pulkovo PUL2 photographic catalog (Bobylev et al. 2004) and HIPPARCOS. 

(5) The KIEV solution. The stellar proper motions from the GPM1 (Rybka and 
Yatsenko 1997) and HIPPARCOS catalogs were compared by Kislyuk et al. (1997). 

(6) The POTSDAM solution. The parameters u x , u> y , u z of the Potsdam program were 
taken from Hirte et al. (1996). 

(7) The BONN solution. The results of the Bonn program are presented in Geffert et 
al. (1997) and Tucholke et al. (1997). 

(8) The EOP solution. The results of the analysis of Earth orientation parameters 
(EOP) were taken from Vondrak et al. (1997). Only two rotation parameters, u x , and uj y 
are determined in this method. 

(9) The HST solution. The results of stellar observations with the Hubble Space 
Telescope (HST) were taken from Hemenway et al. (1997). 

Now, we will point out several new results that have not been used previously to solve 
this problem. 
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(10) The XPM solution. The XPM catalog (Fedorov et al, 2009) contains absolute 
proper motions for about 275 million stars fainter than 12m derived by comparing their 
positions in the 2MASS and USNO-A2.0 (Monet 1998) catalogs. The absolutization 
was made using about 1.5 million galaxies from the 2MASS catalog of extended sources. 
Thus, the XPM catalog is an independent realization of the inertial frame of reference. 
The stellar proper motions from the XPM and UCAC2 (Zacharias et al. 2004) catalogs 
were compared by Bobylev et al. (2010). The parameters u x ,u y ,u z were calculated using 
about 1 million stars. Among all of the programs listed in Table 1, the XPM solution is 
unique in that it was obtained from the differences of stars covering the entire sky almost 
completely, except for the zone 5 > 54°, where there are no UCAC2 stars. 

(11) The VLBI-07 solution. Boboltz et al. (2007) analyzed the positions and proper 
motions of 46 radio stars and obtained new parameters of the mutual orientation of the 
optical realization (HIPPARCOS) and the radio system. 

(12) The "MINOR PLANETS" solution. Chernetenko (2008) estimated the rotation 
parameters of the HIPPARCOS system relative to the DE403 and DE405 coordinate 
systems of ephemerides by analyzing a long-term series of asteroid observations. This 
result suggests that either the dynamical DE403 and DE405 theories need to be improved 
or the HIPPARCOS system needs to be corrected. We reduced the weight of this solution 
by half because of the possible contribution from the inaccuracy of the dynamical DE403 
and DE405 theories. 

The weight of each of the comparison catalogs was taken to be inversely proportional 
to the square of the random error e w in the corresponding quantities u x ,u y ,u z and was 
calculated from the formula 



Not all of the authors use the equations to determine u! x ,ui y , oo z in the form in which they 
were suggested by Lindegren and Kovalevsky (1995): 



where the catalog-HIPPARCOS differences are on the left-hand sides of the equations. 
Therefore, in several cases, the signs of the quoted quantities were reduced to the necessary 
uniform form (Zhu 2001, 2003; Boboltz et al. 2007). 

The last rows of Table 1 give Mean 1 calculated as a simple mean and Mean 2 that 
was calculated as a weighted mean and is the main result of our analysis. 

Denote the components of the rotation vector around the rectangular equatorial axes 
by u x ,u y ,u z ; then, 




(9) 



Afi a cos 5 = u x cos a sin 5 + u y sin a sin 5 — u z cos S, 
Afis = —u x sin a + Uy cos a, 



(10) 
(11) 



( 




a, 



) 



G 



( 



M, 
M] 



21 



32 



13 



) 




(12) 
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where 



G = 



/ -0.0548 
-0.8734 
V -0.4838 



+0.4941 
-0.4448 
+0.7470 



-0.8677* 

-0.1981 

+0.4560 



(13) 



is the well-known transformation matrix between the unit vectors of the Galactic and 
equatorial coordinate systems. From Eqs. (9) and (10), it is easy to see the relationship 

Assuming the components in the mean-2 solution to be zero, 



between M 13 and tu z 




4.74 



I 


\0.52 



(14) 



for the case where the left- and right-hand sides are expressed in km s kpc 



KINEMATICS OF THE DISK WARP 

The results of determining the kinematic parameters of the Ogorodnikov-Milne model 
using a sample of RGC stars derived by simultaneously solving the system of equations 
(6) and (7) are presented in Table 2. The second column gives the solution obtained 
without applying any corrections to the input data; the third column gives the solution 
for the case where the stellar proper motions in the form fx a cos 5 were corrected by 
applying the correction oj z = —0.52 mas yr™ 1 (14) using Eq. (10). 

The parameters derived by separately solving Eqs. (8) and (7) are presented in Table 

3. 

As can be seen from Table 2, applying the correction affected only three components 
of the rotation tensor: insignificantly M~, noticeably M 3 ~, and most strongly M~ 3 , which 
is explained by the structure of the matrix G (13). Since the components M~ and M+ 
found (the third column of Table 2) do not differ significantly from zero, they may be set 
equal to zero. 

Suppose that the values of M 2 ~ and M+ (the Oort constants) found describe only the 
rotation around the Galactic z axis, while the motion in the yz plane is independent. 
Let us now consider the displacement tensor M w that describes the kinematics in the yz 
plane: 



w 



M 22 
M 32 



M 23 
M 33 



/ du 2 du 2 \ 
~dr2~ ~dr^ 
dus <9«3 



(15) 



As has already been noted, when using only the stellar proper motions, we can determine 
only the difference (M+ - M+). The identity (M+ - M+) = (M 33 - M 22 ) is valid for 
the diagonal elements. As can be seen from Table 2, the value of this quantity differs 
significantly from zero. 

Three cases are possible when analyzing tensor (15): 
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1) M 22 ^ 0, M 33 = 0; 

2) M 22 = 0, M 33 ^ 0; 

3) M 22 = 0, M 33 = 0. 

For the completeness of the picture, note that case 4 is also possible: M 22 ^ 0, M 33 7^ 0; 
since this requires data on the stellar radial velocities, this case is not considered here. 

Consider case 1, M 33 = 0,, using the data from the third column of Table 2, M 22 = 
1.3 ± 0.4 km s _1 kpc -1 . The components of the displacement tensor Mw, the symmetric 
deformation tensor Myy, and the antisymmetric rotation tensor are (km s _1 kpc -1 ) 

The deformation tensor in the principal axes is (km s" 1 kpc -1 ): 



Ai 
A 2 
1.8 
-0.5 

and the angle between the positive direction of the Oy axis and the first principal axis of 
this ellipse is 29 ± 6°. 

Equation (1) can be written as 

V = V + grad F + (wxr). 

Analysis of grad F shows (Sedov 1970) that an infinitesimal sphere of radius r composed 
of points of the medium at time t, 

2,2,2 2 
x + y + z = r 

transforms into a deformation ellipsoid after time At, 

x* 2 y* 2 z* 2 9 



(1 + AxAt) 2 (l + A 2 At) 2 (1 + A 3 A*) : 

which is shown in Fig. 2b for our two-dimensional case. 

For case 2, we have M 33 = —1.3 ± 0.4 km s -1 kpc" 1 and M 22 = 0. The deformation 
tensor in the principal axes is (km s -1 kpc -1 ) 



M+ = 



0.5 
-U 
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Table 2: Kinematic parameters of the Ogorodnikov-Milne model found by simultaneously solv- 
ing the system of equations (6) and (7) 



Parameter 


Without correction 


With correction 


x & 


7.99 ±0.10 


7.99 ±0.10 




16.40 ±0.10 


16.40 ±0.10 


z Q 


6.72 ±0.09 


6.72 ±0.09 


A = M+ 


15.82 ±0.21 


15.82 ±0.21 




-1.40 ±0.18 


-2.60 ±0.18 


Mf 3 


-1.99 ±0.18 


-0.15 ±0.18 


B = M 2 1 


-11.99 ±0.15 


-10.87 ±0.15 


Mtl-22 


-7.75 ±0.39 


-7.75 ±0.39 


M+ 


-0.47 ±0.23 


-0.47 ±0.23 




1.00 ±0.22 


1.00 ±0.22 


-^33-22 


-1.26 ±0.44 


-1.26 ±0.44 



Note: Xq,Yq, Z@ in km s , other parameters in km s kpc 



Table 3: Kinematic parameters of the Ogorodnikov-Milne model found by separately solving 
Eqs. (8) and (7) 



Parameter 


Without correction" 


With correction 61 


Without correction b 


With correction 6 


X & 


7.89 ±0.12 


7.89 ±0.12 


8.51 ±0.21 


8.51 ±0.21 


Y Q 


16.13 ±0.13 


16.13 ±0.13 


17.40 ± 0.21 


17.40 ± 0.21 


Z & 






6.73 ± 0.08 


6.73 ± 0.08 


A = M+ 


15.76 ±0.25 


15.76 ±0.25 


15.99 ± 0.48 


15.99 ±0.48 


^32 






-1.47 ±0.25 


-2.66 ±0.25 


M 13 






-1.46 ±0.25 


0.38 ±0.25 


B = 


-12.02 ±0.17 


-10.91 ±0.17 






Kl-22 


-8.21 ±0.46 


-8.21 ±0.46 


-4.31 ± 1.03 


-4.31 ± 1.03 


M+ 






0.17 ±0.34 


0.17 ±0.34 


+ ^2 + 3 






1.11 ±0.32 


1.11 ±0.32 


-^33-22 






0.29 ± 0.59 


0.29 ±0.59 


M 23 


2.24 ±0.46 


3.43 ± 0.46 






M13 


-2.93 ± 0.43 


-1.09 ±0.43 







a The parameters were found by solving Eq. (8). 
fe The parameters were found by solving Eq. (7). 
Note: Xq, Yq, Zq in km s _1 , other parameters in km s _1 kpc 



11 



the angle between the positive direction of the Oy axis and the first principal axis of this 
ellipse is 29 ± 6°. 

For case 3, where M 33 = and M 22 = 0, the deformation tensor has two roots: 
Ai = 1.0 km s -1 kpc -1 and A2 = —1.0 km s _1 kpc -1 ; therefore, the first principal axis of 
this ellipse is oriented at an angle of 45° to the Oy axis. 

In all three cases, the divergence 0.5(M 22 + M 33 ) is insignificant. It is +0.6 ± 0.4, and 
km s" 1 kpc -1 for the first, second, and third cases, respectively. 

As a result, we can conclude that the kinematics in the yz plane can be described as 
a rotation around the Galactic x axis with an angular velocity flw = M~ 2 — X±. This 
angular velocity is —4.4 ± 0.5, —3.1 ± 0.5, and —3.6 ± 0.3 km s -1 kpc -1 for the first, 
second, and third cases, respectively. 

Let us estimate the linear velocities in the yz plane. The mean photometric distance 
for our sample of RGC stars is 7 = 0.57 ± 0.17 kpc. The linear solidbody rotation 
velocity is then M~ 2 ■ r = —1.5 ± 0.1 km s -1 . The maximum deformation velocity is 
Ai • r — 1.0 ± 0.3 km s" 1 and the highest (among the cases considered) linear velocity is 
VL w -r = -2.5 ±0.3 km s _1 . 

Figure 2 gives a qualitative picture that reflects the pattern of our results. The dis- 
tribution of rotation velocity vectors in the yz plane for the case of rotation we found 
is shown in Fig. 2a. Figure 2b shows how a circumference of unit radius turns into a 
deformation ellipse after a unit time interval. 

According to the results of our separate solution of Eqs. (8) and (7) (the third and 
fifth columns of Table 3), for example, for case 1, we have 



We can see from our comparison of matrixes (19) and (16), (20) and (17), (21) and (18) 
that there are differences only in the diagonal element. This leads us to conclude that both 
approaches used (the simultaneous and separate solution of the equations) are equivalent, 
because they yield coincident results. 

When writing Eq. (14), we set the components equal to zero, u x = and oo y = 0. As 
we see from the Mean 2 solution in Table 1, the value of the component u y was found 
to be nonzero outside the 2a range. However, the kinematic parameters calculated by 
taking this into account, i.e., u y 7^ and u z 7^ 0, have no significant differences between 
the solutions reflected in Tables 2 and 3 and, hence, we do not give them. 



Drimmel et al. (2000) provide arguments for the model of precession of the warped disk in 
the zx plane (rotation around the y axis) with an angular velocity of —25 km s -1 kpc -1 . 




(21) 



(19) 



(20) 



DISCUSSION 
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However, as can be seen from Table 2, as a result of applying the correction oo z = —0.52 
mas yr _1 , the component of the rotation tensor around the y axis (M~) has an almost 
zero value and the component of the deformation tensor in this plane (M+) does not differ 
significantly from zero. 

The angular velocity Qw estimated here is inconsistent with the analysis of the motion 
of giants of various spectral types (O-M) performed by Miyamoto et al. (1993), who found 
a positive direction of rotation around the x axis (directed toward the Galactic center). 
Note that the analyzed stellar proper motions were determined in the FK5 (Fundamental 
Catalog 5) system, which is noticeably distorted by the uncertainty in the precession 
constant. On the other hand, having analyzed the proper motions of 0-B5 stars in 
the HIPPARCOS system, Miyamoto and Zhu (1998) also reached the conclusion about a 
positive rotation around the 

At present, it is impossible to choose between the three forms of the deformation tensor 
considered in the previous section due to the absence of data on stellar radial velocities 
and distances. Note that based on a sample of 3632 relatively close RGC stars with known 
HIPPARCOS trigonometric parallaxes and radial velocities, we estimated the parameters 
M 22 = -1.3±3.2 km s" 1 kpc" 1 and M 33 = -1.1 ±2.6 km s" 1 kpc" 1 (Bobylev et al. 2009). 
Consequently, the difference (M+ — M+) is very close to zero and, hence, the orientation 
of the deformation ellipse is close to 45°, as in case 3 considered. 

In the Introduction, we listed various hypotheses about the nature of the Galactic 
disk warp. In our opinion, the hypothesis about an outflow of gas from the Magellanic 
Clouds (Olano 2004) has been worked out in greatest detail. Assuming the distance to 
the Large Magellanic Cloud to be t^mc — 50 kpc, we find that the linear velocity of the 
rotation we revealed at this distance, (—3.1 ± 0.5) < Qw < (—4.4 ± 0.5) km s _1 kpc -1 , 
is \Qw\ ■ f'LMC = (155 4- 220) ± 25 km s _1 . Such velocities are in agreement with the 
estimate of the motion of high- velocity hydrogen clouds relative to the local standard of 
rest, ~ 200 km s" 1 (Olano 2004), and the rotation direction we found is in good agreement 
with the direction of motion of the Magellanic Clouds around the Galaxy during the past 
570 Myr (Fig. 6 from Olano 2004). 

CONCLUSIONS 

We analyzed the kinematics of about 82000 RGC stars from the Tycho-2 catalog. These 
stars are a variety of "standard candles", because they occupy a compact region on the 
Hertzsprung-Russell diagram. Therefore, photometric distance estimates with an accu- 
racy of at least 30% are available for them. The RGC stars considered lie in the range of 
heliocentric distances 0.3 kpc < d < 1 kpc. 

Since these stars are distributed fairly uniformly in space, they are of great interest 
not only in studying the Galaxy's general rotation described by the rotation parameters 
in the Galactic xy plane but also in the other two planes, namely yz and zx. 

For a reliable description of their kinematic peculiarities, we should be confident that 
the observational data (Tycho-2 proper motions) are free from the systematic errors re- 
lated to the referencing of the optical realization of the ICRS/HIPPARCOS system to the 
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inertial frame of reference specified by extragalactic sources. 

Here, we gave much attention to this problem. Based on all of the currently available 
data, we determined new, most probable components of the residual rotation vector of 
the optical realization of the ICRS/HIPPARCOS system relative to the inertial frame of 
reference, (u x ,u y ,u z ) = (-0.11,0.24,-0.52) ± (0.14,0.10,0.16) mas yr -1 . This led us to 
conclude that a small correction, u z = —0.52 mas yr -1 , should be applied to the Tycho-2 
stellar proper motions of the form /j, a cos 5. 

By applying this correction, we showed that, apart from their involvement in the 
general Galactic rotation described by the Oort constants A = 15.8 ± 0.2 km s -1 kpc -1 
and B = —10.9 ± 0.2 km s -1 kpc -1 , the RGC stars considered have peculiarities in the 
yz plane related, in our opinion, to the kinematics of the warped Galactic stellar-gaseous 
disk. 

The component of the solid-body rotation of the local solar neighborhood around the 
Galactic x axis is M~ 2 = —2.6 ± 0.2 km s _1 kpc -1 . The parameter of the deformation 
tensor in this plane M+ = 1.0±0.2 km s -1 kpc -1 and the difference M 33 -M 22 = -1.3±0.4 
km s -1 kpc -1 differ significantly from zero. 

On the whole, the kinematics of the warped Galactic stellar-gaseous disk in the local 
solar neighborhood can be described as a rotation around the Galactic x axis directed 
from the Sun toward the Galactic center with an angular velocity (—3.1 ± 0.5) < Qw < 
(—4.4 ±0.5) km s -1 kpc -1 . Thus, the rotation is around an axis close to the line of nodes 
of this structure. In the case where the difference M 33 — M 22 may be set equal to zero, 
the angular velocity Vt w is —3.6 ± 0.3 km s -1 kpc -1 . If M 22 = and M 33 ^ 0, then 
VL W = -3.1 ± 0.5 km s -1 kpc -1 . If M 22 ^ and M 33 = 0, then Q w reaches -4.4 ± 0.5 km 
s -1 kpc -1 . 
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